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Abstract 



We report on the correlation between the flux, color and polarization variations on time scales of days- 
months in blazars, and discuss their universal aspects. We performed monitoring of 42 blazars in the optical 
i-C ' and near-infrared bands from 2008 to 2010 using TRISPEC attached to the "Kanata" 1.5-m telescope. 

We found that 28 blazars exhibited "bluer-when-brighter" trends in their whole or a part of time-series 
data sets. This corresponds to 88% of objects that were observed for > 10 days. Thus, our observation 
unambiguously confirmed that the "bluer-when-brighter" trend is common in the emission from blazar 
jets. This trend was apparently generated by a variation component with a constant and relatively blue 
color and an underlying red component. Prominent short-term flares on time scales of days-weeks tended 
to exhibit a spectral hysteresis; their rising phases were bluer than their decay phases around the flare 
Cn| \ maxima. In contrast to the strong flux-color correlation, the correlation of the flux and polarization degree 

was relatively weak; only 10 objects showed significant positive correlations. Rotations of polarization were 
detected only in three objects: PKS 1510-089, 3C 454.3, and PKS 1749+096, and possibly in S5 0716+714. 
We also investigated the dependence of the degree of variability on the luminosity and the synchrotron 
peak frequency, f pC ak- As a result, we found that lower luminosity and higher i/ pca k objects had smaller 
variations in their amplitudes both in the flux, color, and polarization degree. Our observation suggests 
the presence of several distinct emitting sources, which have different variation time-scales, colors, and 
polarizations. We propose that the energy injection by, for example, internal shocks in relativistic shells is 
a major factor for blazar variations on time scales of both days and months. 

Key words: galaxies:active — galaxies: BL Lacertae objects: general — galaxies: jets 



. 1. Introduction tering of synchrotron emission and/or external photons 

(e.g. Dermer et al. 1992; Sikora et al. 1994). 
Blazars are a subgroup of active galactic nuclei (AGN). Blazars are classified into two categories: BL Lac ob- 
They have relativistic jets that are believed to be directed jects and flat spectrum radio quasars (FSRQs) (Angel, 
along the line of sight (e.g. Blandford, Rees 1978). The Stockman 1980). Non-thermal continuum emission is 
emission from blazars can be detected in a very wide range dominant in BL Lac objects, while broad emission lines 
of wavelengths from radio to TcV 7-ray regions, with from AGN characterize FSRQs. BL Lac objects are fur- 
strong variability on various time-scales (e.g. Hufnagel, ther classified into three subtypes based on the peak frc- 
Bregman 1992). Their spectral energy distribution (SED) quency of synchrotron emission: high-, intermediate-, and 
is characterized by non-thermal continuum spectra that low-energy peaked BL Lac objects (HBL, IBL, and LBL). 
consist of low and high-energy components. The low- The peak frequency lies higher than, around, and lower 
energy component is believed to be synchrotron radia- than the optical band in HBL, IBL, and LBL, respectively, 
tion from relativistic electrons in the jet. The emission Multi-wavelength observations are essential to study 
of this component is highly polarized (Angel, Stockman blazars because of their emission in a wide range of wave- 
1980; Mead et al. 1990). The origin of the high-energy lengths. The multi- wavelength study of blazars has re- 
component is not fully understood. The most plausible cently entered a new era owing to 7-ray observations ob- 
scenario is that it is emission via inverse Compton scat- tained with the Fermi satellite (Abdo et al. 2010a; Abdo 
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et al. 2010b). The origin of the high-energy component, 
for example, is expected to be revealed by the study of 
temporal variations in SEDs obtained with Fermi. Among 
the multi-wavelength study, optical data is important to 
understand the SED variations, because the optically-thin 
synchrotron emission in this band allows us to estimate 
basic physical parameters in the emitting region. For 
follow-up ground-based observations of Fermi blazars, we 
executed photopolarimetric monitoring of 42 blazars from 
2008 to 2010. We obtained simultaneous optical and near- 
infrared (NIR) data using the "Kanata" 1.5-m telescope. 
Our observation provides one of the largest data sets of 
blazars in terms of variations in color and polarization on a 
time scale of days-months (Ikejiri et al. 2009). In this pa- 
per, we report on the results of simple correlation studies 
between the flux (,or luminosity), color, polarization de- 
gree (PD), and polarization angle (PA). We also discuss 
the implications for the mechanism of time variations in 
blazars obtained from our observations. A detailed study 
of polarization and the data, itself, will be published in a 
forthcoming paper. 

The behavior of colors gives us a clue to understand 
the mechanism of time variations in blazars. Carini et al. 
(1992) report a possible feature that both BL Lac and 
OJ 287 became bluer when they were brighter (also see, 
Clements, Carini 2001). Ghiscllini et al. (1997) observed 
S5 0716+714, and found that a bluer-when-brighter trend 
was seen only in its faint state. Its bluer-when-brighter 
trend was later confirmed in variations on a time scale 
shorter than days in its bright state (Wu et al. 2007a; 
Sasada et al. 2008). Villata et al. (2002) and Villata 
et al. (2004) report that the bluer-when-brighter trend 
in BL Lac was only observed in short-term variations, 
and not prominent in long-term ones. Raitcri et al. 
(2001) performed multi-band photometric observations of 
AO 0235+164 for four years, and found a bluer-when- 
brighter trend. As reported in those past studies, the 
bluer-when-brighter trend has been a well-observed fea- 
ture in blazars, while its universality has not been estab- 
lished. 

Systematic observations of multiple sources have also 
been performed in order to establish a characteristic fea- 
ture of color variations in blazars. Ghosh et al. (2000) 
performed observations of five blazars. Their observation 
showed that 3C 66A only exhibited a bluer-when-brighter 
trend, while the other objects showed no significant cor- 
relation of the flux and color. Gu et al. (2006) inves- 
tigated variations in color of eight blazars. The blucr- 
when-brighter trend was confirmed in five blazars. Among 
the other three objects, 3C 345 showed no correlation 
of the flux and color, and 3C 454.3 and PKS 0420-01 
showed a reddening trend when they brightened. This 
"redder-when-brighter" trend in 3C 454.3 was confirmed 
in later observations (Raiteri et al. 2008b, Villata et al. 
2006a). Thus, it is currently unclear whether the blucr- 
when-brighter trend is universal in blazars. 

As well as the color variation, the polarimetric vari- 
ation in the optical band has been extensively investi- 
gated in blazars. The temporal variation in polarization 



has, in general, been considered to be erratic in blazars 
(e.g. Moore et al. 1982). On the other hand, a few cases 
have been reported in which the increase in PD was as- 
sociated with flares of the total flux. For example, PD of 
Mrk 421 increased to ~ 14 % associated with its outburst 
in 1997 (Tosti et al. 1998). A significant correlation of 
the flux and PD was also detected in AO 0235+164 in 
2006 (Hagen-Thorn et al. 2008). Most recently, several 
rotation episodes of PA have been found during flares 
(Marscher et al. 2008; Marschcr ct al. 2010; Abdo et al. 
2010c; Jorstad et al. 2010). In general, the number of po- 
larimetric observations has been much smaller than that 
of color observations. High-frequency, long-term polari- 
metric observations are required in order to find universal 
aspects that are possibly present. 

The aim of this paper is to investigate universal as- 
pects in blazar variability in the optical-NIR bands. We 
search for them with a simple correlation study between 
the flux, color, and polarization. In section 2, we de- 
scribe the observation method and reduction processes. 
In section 3, first, we shortly introduce the basic proper- 
ties of our sample (subsection 3.1). Then, we report on 
the results of our correlation studies of the flux and color 
(subsection 3.2), and the flux and polarization (subsec- 
tion 3.3). We also search for rotation events of polariza- 
tion in subsection 3.3. In subsection 3.4, we report on 
the dependence of the degree of variability on different 
luminosities and synchrotron peak frequencies. In sec- 
tion 4, we discuss the implication to the mechanism of 
blazar variability obtained from our observation. Finally, 
we summarize our findings in section 5. 

2. Sample and Observation 

2. 1 . Sample 

We first selected our targets from the catalog, 
"Extended list of 206 possible AGN/blazar targets for 
GLAST multi-frequency analysis" 1 with an apparent 
magnitude of R < 16. Then, we performed test obser- 
vations of those potential targets in early 2008. Objects 
were included in our sample in the case that the test obser- 
vation confirmed that they were bright enough to monitor 
with the photopolarimetric mode (R < 16). The number 
of selected objects from that catalog was 30. New objects 
were included in our sample when optical or 7-ray flares 
of them were detected during our monitoring period. The 
number of additional objects was 12. The total number of 
our target is, hence, 42. Table 1 lists our targets. Bright 
and/or highly variable sources were monitored more fre- 
quently than faint and/or stationary ones. This is because 
one of our aims was follow-up observations of 7-ray flares 
detected by Fermi. 

The objects were classified into four subclasses, de- 
fined in Abdo et al. (2010a). Based on high-quality 
SEDs obtained almost simultaneously from radio to 7- 
rays, Abdo et al. (2010a) propose a new classification for 
blazars (FSRQs and BL Lac objects): Low, intermedi- 

1 http:/ /glastweb. pg.infn.it/blazar/ 
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Table 1. Observation log and list of comparison stars. 
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(l)Object name. (2)Period of our monitoring. (3)Number of observations. 

(4)Exposure times for each V, J, and iC s -band image. (5)Coordinate of comparison stars. 

(6),(7),(8)V , J, and Ks band magnitudes of comparison stars. 

(9) Rcfcrences for the magnitudes of comparison stars. 

(10) Standard deviations of the V-mag of the comparison stars during our monitoring period. The magnitudes were measured 
with neighbor check stars in the same field. Note that the check stars are fainter than the blazars in some cases, 

and hence, cry do not indicate the systematic error of the photometry of the blazars. It just 
indicates the constancy level of the comparison stars. 

[l]Skiff (2007); [2]Villata et al. (1998); [3] Gonzalez- Perez et al. (2001); [4] Adelman-McCarthy et al. (2008); 

[5]Doroshenko et al. (2005); [6]H0g et al. (2000); [7]McGimsey, Miller (1977). 

*We measured the V-magnitude of our comparison star using the stars listed in the reference. 

fThc NIR magnitudes were referred from the same reference for the V-band magnitude. 

No symbol: The V-magnitude of the comparison star was quoted from the reference, and the NIR 

magnitudes were from the 2MASS catalog. 



ate, and high synchrotron peaked blazars (LSP, ISP, and 
HSP blazars, respectively). LSP, ISP, and HSP blazars 
are defined as those having a synchrotron peak frequency, 
^pcak, of i/ poak < 10 14 Hz, 10 14 < z/peak < 10 15 Hz, and 
^peak ^ 10 15 Hz, respectively. All FSRQs are LSP blazars 
in their sample, except for one source (J1012. 9+2435; 
ISP), which is not included in our sample. BL Lac ob- 



jects are distributed in all three subclasses. According to 
Abdo et al. (2010a), the mean ^ poa k of LSP BL Lac objects 
is higher than that of FSRQs. This indicates that ^ pe ak 
is lower in the following order: FSRQs < LSP-BL Lac < 
ISP-BL Lac < HSP-BL Lac. We classified blazars into 
these four subclasses, namely FSRQs, LSP BL Lac, ISP 
BL Lac, and HSP BL Lac, and used abbreviated forms of 
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Table 2. Observational properties of blazars. 
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the latter three subclasses as LBL, IBL, and HBL. 

We identified 27 of our monitoring objects in the cata- 
log presented in Abdo ct al. (2010a). The classification of 
these 27 objects was based on Abdo et al. (2010a). Among 
the remaining 15 objects, prominent emission lines have 
been observed in QSO J0324+3410 (Marcha et al. 1996; 
Zhou et al. 2007), PKS 1222+216 (Burbidge, Kinman 
1966), 4C 14.23 (Healey et al. 2008), and 3EG J1236+0457 
(Halpcrn ct al. 2003). Hence, these four objects belong 



to FSRQs. Then, the classification of the remaining 11 
objects was based on Nicppola et al. (2006). According 
to Nicppola ct al. (2006), BL Lac objects are classified 
into three subgroups: namely LBL, IBL, and HBL, de- 
fined by t'peak estimated from radio-optical (with X-rays 
in several cases) SEDs. Their criteria for subgroups are, 
however, different from those in Abdo ct al. (2010a): 
tVak - 1014 5 Hz for LBL ' 1014 5 - ^pcak < iO 16 ' 5 Hz for 
IBL, and is pca , k > 10 16 5 Hz for HBL. We classified the 



No. ] 

remaining 11 objects by adopting the criteria defined in 
Abdo et al. (2010a) using ^ pca k estimated in Nieppola 
et al. (2006). Table 2 gives v pea k and subclasses of each 
object. Our sample includes 13 FSRQs, 8 LBLs, 9 IBLs, 
and 12 HBLs. We note that ^ pca k can significantly shift, 
depending on the brightness of blazars (e.g. Costamante 
et al. 2001). It is possible that the classification of each 
object changes with time, even with the same criterions. 
Table 2 also presents the red-shift, z, and the variation 
range obtained from previous literature. 

2.2. Observation and data reduction 

We performed photopolarimctric observations of those 
blazars simultaneously in the optical V-, and the NIR J-, 
and if s -bands from 2008. We used TRISPEC attached 
to the 1.5-m "Kanata" telescope at Higashi-Hiroshima 
Observatory. TRISPEC is capable of simultaneous three- 
band (one optical and two NIR bands) imaging or spec- 
troscopy, with or without polarimctry (Watanabc et al. 
2005). TRISPEC has a CCD and two InSb arrays. Table 1 
presents the observation period, the number of observa- 
tions, and the typical exposure times of each image. A 
part of NIR images failed to be obtained because of me- 
chanical errors. The if s -band data are, especially, limited, 
and have a low quality. In this paper, we mainly use V- 
and J-band data for the color and correlation analysis, 
except for in sub-subsection 3.2.1. 

The photometry of blazars was performed in a standard 
procedure of CCD images; after making dark-subtracted 
and flat-fielded images, the magnitudes were measured us- 
ing the aperture photometry technique. The radius of the 
aperture, which depended on the seeing size of each night, 
was 3-5 arcsec. These correspond to 3-4 pixels on the op- 
tical CCD. Since the pixel scale of the NIR InSb arrays of 
TRISPEC is large (~ 1.66 arcsec), the NIR images were 
slightly de-focused in order to avoid undersampling. On 
the NIR arrays, the FWHM of an point source was ad- 
justed to be ~ 3 pixels, which correspond to the aperture 
size of the NIR images. It is well known that the photom- 
etry of a part of blazars is complicated by the presence 
of a host galaxy component. Subtraction of contamina- 
tion of the host galaxy is, in general, possible if the host 
galaxy is well resolved on images (e.g. Kotilainen et al. 
1998; Nilsson et al. 2007). Unfortunately, the angular res- 
olution of our instrument is too low to accurately estimate 
the host galaxy contamination. Hence, we performed no 
correction for it in the present work. The obtained magni- 
tude corresponds to the combined flux from the jet, AGN, 
and the host galaxy. 

We calculated differential magnitudes of blazars using 
comparison stars located in the same frame. We checked 
the constancy of the brightness of the comparison stars us- 
ing the differential photometry between them and neigh- 
bor stars in the same field. The standard deviations of 
the V-band magnitudes, cry, are given in table 1. We 
confirmed that each comparison star exhibited no signifi- 
cant variation. The V-band magnitude of the comparison 
stars were taken from previous literature. The references 
are given in the ninth column of table 1. The symbol "*" 



5 

represents cases that our image includes no comparison 
stars listed in the reference. In this case, we measured 
the V-magnitudc of our comparison star using the stars 
listed in the reference. The J- and i^ s -band magnitudes 
were referred from the 2MASS catalog (Skrutskie et al. 
2006), except for the cases with the symbol "f" in table 1. 
The NIR magnitudes with this symbol were taken from 
the same reference for the V-band magnitude. 

The correction for the interstellar extinction, Ay, and 
the color excess were performed according to Schlegcl et al. 
(1998). The magnitude and color presented in this paper 
are corrected values. The correction was not done for the 
historical variation range given in table 2. 

A set of polarization parameters was calculated from 
four consecutive images, which were obtained with half- 
wave-plate angles of 0°, 22.5°, 45.0°, and 67.5°. We 
took 12 sets of images for each object on one night, from 
which three sets of polarimctric data were obtained. We 
confirmed that the instrumental polarization was smaller 
than 0.1 % in the V-band using observations of unpolar- 
ized standard stars. Hence, we applied no correction for 
it. The zero point of the polarization angle was corrected 
as the standard system (measured from north to east). 

Since we present no discussion about intra-day variabil- 
ity in this paper, we only use nightly-averaged photomet- 
ric and polarimctric data. Observations were sometimes 
carried out under bad sky conditions. Some data obtained 
under such conditions have very large errors. They could 
just disrupt systematic trends, which may exist in blazar 
variability. In this paper, we only use photometric data 
with an error of less than 0.1 mag, color indices with an 
error less than 0.1, and PD with an error less than 5%. 
In Appendix, we show the temporal variation in the flux, 
color, PD, PA, and the color-magnitude and magnitudc- 
PD diagrams, and the Stokes QU plane for all objects. 
Those data will be published in electric form in a forth- 
coming paper about a detailed analysis of polarization 
(Paper II). 

3. Results 

3.1. Basic Properties of the Sample 

3.1.1. Relationship between the red-shift, luminosity, and 
synchrotron peak frequency, with comments on the 
blazar sequence 

In this subsection, we report on basic properties of our 
sample. We investigated the relationships between the 
red-shift, z, the optical luminosity, vL v , and the syn- 
chrotron peak frequency, z/ pC ak, of our sample. We cal- 
culated the minimum, average, and maximum vh v dur- 
ing our observation period for our sample whose z is 
known. The observed minimum and maximum values 
of the V-band magnitudes and z are presented in ta- 
ble 2. No estimation of vL u was given in the case that 
z has not been well determined. Such z are indicated in 
parentheses in table 2. We assumed a flat universe with 
Hq = 71 kms -1 Mpc -1 and fl m — 0.27 for calculating vL u . 

Figure 1 shows the z distribution of our sample. Our 
sample is distributed in z = 0.02-1.75. The vL v rep- 



Photopolarimetric Monitoring of Blazars 



G 



Y. Ikejiri et al. 




^10 44 



v 10" 




10 



,43 



-32 
-30 
-28 

-26 > 
-24 
-22 
-20 



_l I I L: 



12 13 14 15 16 17 18 19 20 

l0 9 v peak ( Hz ) 



Fig. 1. z distribution of the V-band luminosity (vL v ) (and 
absolute magnitude, My) of our sample. The circles and tri- 
angles represent the average vL v during our observation pe- 
riod for BL Lac objects and FSRQs, respectively. The filled 
and open circles represent the observed and A"-corrected lumi- 
nosities, respectively. The vertical bars associated with each 
point show variation ranges defined by the observed minimum 
and maximum values of the magnitude during our observation 
period. The dashed line indicates an apparent magnitude of 
y = 17.0 

resented by the open circles were calculated from the 
observed magnitudes with a A"-correction. For the K- 
correction, we assumed SEDs in two power-law forms be- 
tween the V — J-band and J — A s -band regions. The flux 
densities at 0.55, 1.22, and 2.16 /im were obtained from 
our simultaneous V-, J-, and A s -band observations, re- 
spectively. Then, the rest-frame flux densities at 0.55 /im 
were estimated with those three flux densities and the two- 
powcr-law SEDs. The dashed line in figure 1 represents 
an apparent magnitude of V = 17.0. This curve indicates 
that our sample had a limit magnitude of V ~ 17. 

Figure 2 shows the ^ poa k distribution of the luminos- 
ity. In addition to the AT-corrcction, a red-shift correction 
was performed to f pea k of the open circles by multiplying 
k'peak by (1 + z). An anti-correlation is apparent between 
vL v and ^ pC ak- The anti-correlation is more emphasized 
between the corrected vL v and ^ pC ak, as indicated by the 
open circles. 

The anti-correlation between vL u and f pca k is well 
known as the "blazar sequence" , which was first proposed 
by Fossati et al. (1998). According to Ghiscllini et al. 
(1998), the blazar sequence can be explained by the idea 
that a synchrotron cooling works more severely for low- 
energy electrons in more powerful blazars. As a result, 
it leads to smaller f pea k in them. However, it has been 
suspected that the blazar sequence scenario involves se- 
rious problems about the selection effects (Costamante 
et al. 2001; Padovani et al. 2003; Caccianiga, Marcha 
2004; Anton, Browne 2005; Nieppola et al. 2006; Padovani 
2007). When the blazar sequence is evaluated, a prob- 
lem is the high variability not only in vL v , but also 
J^eak (e.g. Costamante et al. 2001). Simultaneous mea- 
surements of vL v and ^ poa k are, hence, required. The 



Fig. 2. i^p Ca k distribution of the V-band luminosity {vL^) 
(and absolute magnitude, My) of our sample. The circles 
and triangles represent BL Lac objects and FSRQs, respec- 
tively. The filled and open circles represent the observed and 
the J£T-corrected luminosities, respectively. The redshift-cor- 
rection was also performed for ^ pea k of the open circles. The 
vertical bars associated with each point show the minimum 
and maximum values of the magnitude during our observa- 
tion period. The dashed line represents the luminosity of a 
typical host galaxy (Urry et al. 2000). The dotted line rep- 
resents an apparent magnitude of V = 17 obtained from the 
z— z/peal; relation in figure 3. 

^ P cak of our sample are mostly quoted from Abdo et al. 
(2010b), who analyzed multi- wavelength data taken in the 
almost same observation periods as those of our monitor- 
ing. Hence, our data can provide one of the most reliable 
tests for the blazar sequence scenario in terms of simulta- 
neous measurements of vL v and J/ pea k- 

As shown in figure 2, the result is apparently compatible 
with the blazar sequence scenario originally proposed by 
Fossati et al. (1998). The blazars with a high z^ poa k of 
> 10 16 Hz had My ~ 6 mag fainter than those with a low 
^ P oak of < 10 Hz. It is evident that our observation did 
not cover the faintest states in several blazars, as can be 
seen in table 2. The observed faintest states in a part 
of LBLs and FSRQs are, in particular, ~ 3 mag brighter 
than the previously recorded faintest states. However, the 
anticorrclation between vL v and f pe ak still appears even if 
vL v are estimated from the faintest states ever recorded. 
This indicates that the flux variability actually has no 
major affect on the blazar sequence scenario. 

We comment on the selection effects. Several studies 
have reported that deep radio observations found sources 
that have low vL v and low j/ pea k (Padovani et al. 2003; 
Caccianiga, Marcha 2004; Anton, Browne 2005). Similar 
results were also obtained in the optical data collected by 
Nieppola et al. (2006). The presence of sources having 
low vL v and low f pca k contradicts the blazar sequence 
scenario. Figure 3 shows the ^ poa k distribution of z of our 
sample. IBLs and HBLs having f pea k > 10 14 Hz can be 
found only in a low z region of z < 0.5. In contrast, a part 
of LBLs and FSRQs can be found in a high z region of 
z > 0.5, as well as in a low z region. The dashed line in 
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Fig. 3. !^p Ca k distribution of the red-shift, z of our sam- 
ple. The circles and triangles represent BL Lac objects and 
FSRQs, respectively. The dashed line represents the best-fit- 
ted power-law model of z against i^p ea k. 

figure 3 represents the best-fitted power-law model of z 
against v-peak- While the power-law form has no physical 
meaning here, this empirical model well represents a trend 
of z versus i/poak- In conjunction with the curve, which in- 
dicates the apparent magnitude of V = 17 in figure 1, this 
model for z provides vL v corresponding to V — 17 in fig- 
ure 2. The dotted line in figure 2 represents it. Objects are 
detectable over this curve in a magnitude-limited sample. 
As suggested by this curve, the apparent lack of low vL v 
and low f p0 ak sources could be due to a selection effect 
present in magnitude-limited samples. 

Second, the lack of high vL v and high z^ poa k blazars is 
also suspected to be due to a selection effect. Giommi 
et al. (2005) propose that such blazars, if exist, would 
have featureless optical spectra, which make it difficult 
to estimate z and luminosity (also sec, Giommi ct al. 
2002). Urry et al. (2000) have reported that host galax- 
ies of BL Lac objects have an absolute magnitude of 
M R = -23.7 ± 0.6. Assuming V - R = 0.7 for elliptical 
galaxies, we show the 1^-band absolute magnitude of the 
typical host galaxy in figure 2 by the dashed line. It 
has been reported that the luminosity of the host galaxy 
is not severely dependent on z in a region of z < 0.5 
(Scarpa et al. 2000; Urry et al. 2000). Our sample in- 
cludes only a few IBLs and HBLs having luminosities one 
order higher than that of the typical host galaxy. A signif- 
icant part of BL Lac objects has such high luminosities in 
"peak ^ 10 14 Hz. They are, however, LBLs in z > 0.5, where 
the absolute magnitudes of the host galaxy are poorly 
known. Thus, our observation is compatible with the sce- 
nario proposed by Giommi ct al. (2005). 
3.1.2. Relationship between the color and synchrotron 
peak frequency 

Figure 4 shows the ^ pca k distribution of V — J. The 
spectral index between the V- and J-band region, otyj, 
defined with f v cx v~ avj , is also shown in the right side 
of the figure. We calculated avj from the flux density at 
0.55 and 1.22 fim estimated from the V- and J-band obser- 
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Fig. 4. ^p C ak distribution of the V — J color index (, and 
the spectral index, ctyj) of our sample. The filled and open 
circles represent the observed and the red-shift corrected data, 
respectively. The vertical bars associated with each point 
show the minimum and maximum values of V — J during our 
observation period. 

vations, assuming a power-law form of spectra. The figure 
indicates an anti-correlation between V — J and f pea k in 
a region of i^pcak ^ 10 16 Hz. This anti-correlation can be 
naturally expected if the optical-NIR emission is dom- 
inated by synchrotron emission from jets. The overall 
trend of V — J (, or ayj) is qualitatively analogous to 
that of the spectral index between the radio and optical 
regions (Fossati et al. 1998; Nicppola et al. 2006). 

In the lowest f pca k regime, objects are distributed 
around oty j = 1.5. This is consistent with the model 
of synchrotron self-Compton emission in blazar jets; the 
spectral index is predicted to be 1.5 in a high energy re- 
gion where Compton cooling is efficient (Chiang, Bottcher 
2002). The optical band definitely corresponds to this 
high energy region in low ^ pca k objects. On the other 
hand, most objects can reach ayj > 1.5, as can be seen 
from the figure. This might suggest that the Compton 
cooling actually worked more efficiently than the model 
predicted. Alternatively, it may suggest that the num- 
ber of high-energy electrons is much smaller than that 
expected from the standard power-law form of the energy 
distribution of electrons. 

The anti-correlation of V — J and f pca k appears to 
be weak in a high i^ pca k region of f pca k ^ 10 15 Hz. In 
this f pea k region, the spectral index is distributed around 
~ 1.0. There is no object in our sample in which fully 
self-absorbed synchrotron emission (a ^ — 1) is domi- 
nant in the optical-NIR band. These colors of high 
Vpeak sources imply that blazars have nearly flat SEDs 
(a ~ 1) in a region even 4-5 orders of magnitude below 
^pc a k- Alternatively, the contamination of the red host 
galaxy may be significant, particularly in low-luminosity 
HBLs. In fact, we found that two low luminosity HBLs, 
1ES 2344+514 (^ pcak = 10 16 4 Hz) and Mrk 501 (i/ pca k = 
10 17 - 1 Hz), had quite red colors of V - J ~ 2.8 and - 2.5, 
respectively, which suggest a large contribution of their 
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host galaxies in the observed fluxes. In addition, it has 
been reported that the UV-IR emission of Mrk 501 can 
be reproduced with a strong contribution of an elliptical 
galaxy (Katarzyhski et al. 2001). 

Before closing this subsection, we mention the anoma- 
lies in a low ^ pca k region in figure 4. First, the red- 
dest source is AO 0235+164 (V - J = 3.1). This ob- 
ject is suggested to be affected by absorption by fore- 
ground galaxies (e.g. Yanny et al. 1989). The quite red 
color of AO 0235+164 is, hence, due to color excess aris- 
ing from the foreground galaxies. According to Raiteri 
et al. (2005), the color excess, E(V — J), arising from 
them is estimated to be 0.673. The mean corrected V — J 
is, then, calculated as 2.44, which is rather normal for 
FSRQs and LBLs having ^ poak - 10 13 ~ 14 Hz. The sec- 
ond reddest source is MisV1436 (V — J = 2.9), whose vi- 
olent variability has been first reported in 2008 by the 
MISAO project. 2 No detailed study has been performed 
for this object. Its atypically red color suggests that this 
blazar is a noteworthy object. Third, the bluest source 
is 3C 273 (V — J = 1.1), one of the most famous QSOs. 
A large contamination of the big blue bump originated in 
its AGN component has been confirmed in optical spec- 
troscopy and polarimetric observations (Impey et al. 1989; 
Smith et al. 1993). The blue color of 3C 273 is due to the 
strong contribution of the AGN component. 

3.2. Light Curve and Color Variation 

3.2.1. Time lag between variations in the optical and 
NIR bands 

Multi-wavelength studies of blazar variability have 
shown that variations in different wave-bands occasion- 
ally correlate with significant time-lags (e.g. Ulrich et al. 
1997). It has been reported that the variation in low- 
energy synchrotron photons tends to lag behind that in 
high-energy ones. In HBLs, for example, variations in 
soft X-rays lag those of hard X-rays by hours (Takahashi 
et al. 1996; Kataoka et al. 2000) and variations in UV 
photons lag those of X-rays or 7-rays by 1-2 d (Buckley 
et al. 1996; Urry et al. 1997). A time lag of low-energy 
photons indicates a spectral hysteresis in flares; for exam- 
ple, a spectrum becomes hardest before the flare maxi- 
mum (Kataoka et al. 2000). Such spectral hystereses have 
also been reported in optical flares (Ciprini et al. 2003; 
Wu et al. 2007b). The time lag and spectral hysteresis 
are expected in the case that synchrotron cooling governs 
temporal variations in a flare region. High-energy elec- 
trons diminished more quickly than low-energy ones by 
cooling (Takahashi et al. 1996; Kirk et al. 1998). 

We searched for universal time lags between the opti- 
cal and NIR variations for each object. We calculated the 
correlation function between the V-, J- and -Kg-band light 
curves. We used the discrete correlation function (DCF) 
method (Edelson, Krolik 1988). The DCF method was 
proposed to estimate a correlation function from unevenly 
sampled data. For two data streams, and bi, the DCF 
method first calculates a set of unbinncd discrete corrcla- 
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Fig. 5. DCF between the light curves in the V- and J-bands 
(panel a), the V- and Ks-bands (panel b), the V-band light 
curve and the V — J color (panel c), and the V-band light 
curve and the V — K s color (panel d) calculated from the 
data of BL Lac. 
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where a and b are the averages of the data values, o~ a and 
(7b are their standard deviations, and e a and ef, are the 
observation error of each data. DCF for each time lag, 
r, is defined as an average of UDCF having the same r: 
Averaging over the M pairs for which r — At/2 < At^ < 
r + At/2, 

DCF(t) = ±J2 UDCF v- ( 2 ) 
The standard error for each bin is: 

°{t) = {YftDCFy - DCF{r)YY/\ (3) 

The upper panels in figure 5 show the DCF between 
the light curves in the V- and J-bands (panel a) and the 
V- and i5T s -bands (panel b) calculated from the data of 
BL Lac. We also calculated the DCF between the V- 
band light curve and the V — J or V — K s color, as shown 
in the lower panels in figure 5. We found no significant 
time-lag in these results for BL Lac. As well as BL Lac, we 
calculated and checked DCF for all blazars we observed. 
We confirmed that no object exhibited a universal signif- 
icant time-lag between variations in the different bands 
and the colors. 

Then, we searched for time lags in the optical and NIR 
light curves for each prominent flare, while no clear lag 
was found between V- and J-band light curves of flares 
in all objects. On the other hand, we found several cases 
that the bluest phase in V — J preceded the flare max- 
imum. Panels (a) and (b) of figure 6 present examples 
for S5 0716+714 and PKS 0048-097, respectively. The 
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Fig. 6. Examples of flares showing time lag between different 
wave-bands. Panel a, b, and c: The V-band light curves 
(filled circles) and V — J color (open circles) during flares in 
S5 0716+714, PKS 0048-097, and S2 0109+224, respectively. 

lags were 2 and 1 d for S5 0716+714 and PKS 0048-097, 
respectively. We systematically searched for time lags be- 
tween the flux and color during such flares. It is expected 
that a time lag possibly present in small flares would read- 
ily disappear by the composition of another small flares. 
Hence, prominent flares should be selected. We selected 
such flares with the following criteria: (i) A flare maxi- 
mum was defined as the brightest point in a period of 15 d 
before and after the maximum date, (ii) A peak-to-peak 
amplitude of the flare must be larger than V = 1.0 within 
30 d. (iii) A flare maximum must be recorded within 5 d 
after and before the neighboring observations. The selec- 
tion with these criteria yielded 17 flares. Among those 
17 events, the peak of the flux coincided with the bluest 
time in nine events. Seven events showed the precedence 
of the bluest phase against the flare maximum. One event 
in S2 0109+224 possibly showed a lag of the bluest time 
behind the flare maximum, which is shown in panel (d) 
of figure 6. All lags were 1-2 d. These results suggest 
that the bluest phase of flares tended to precede the flare 
maximum by < 1 d. 



The spectral hysteresis associated with flares is again 
discussed in sub-subsection 3.2.3. In the next subsection, 
we investigate the correlation between the flux and color, 
neglecting time lags between them. 

3.2.2. Correlation between the V -band light curve and 
the V — J color index 
We present examples of the light curves, color varia- 
tions, and color-magnitude diagrams in figure 7. Panel (a) 
shows those of 3C 371. The color- magnitude diagram 
clearly indicates that the object became bluer when it was 
brighter. In this paper, we conclude that an object had 
a blucr-when-brighter trend in the case that it showed a 
significant positive correlation between the F-band mag- 
nitude and the V — J color index. We used the Pearson 
product-moment correlation coefficient. The significance 
of the correlation was tested using the Student's f-test. 
The correlation coefficient and its 95 % confidence inter- 
val of 3C 371, for example, were calculated as 0.92l ; 04 . 
Thus, we conclude that 3C 371 had a bluer-when-brighter 
trend. We calculated the correlation coefficients using the 
whole data sets of 42 blazars, as well as 3C 371. The 
results are listed as r co i or in table 2. We found that 23 
blazars had significant bluer-when-brighter trends. The 
correlation coefficients that are not statistically significant 
are shown with parentheses in table 2. 

Two blazars, 3C 454.3 and PKS 1502+106, had sig- 
nificant negative correlations between the light curve 
and color, in other words, a redder-when-brighter trend 



(rcoior = -0.67±8;jj?and 



-0.35_ ; 21 , respectively, as shown 
in table 2). Panel (b) of figure 7 shows the light curve 
and color-magnitude diagram of 3C 454.3. As reported 
in Sasada et al. (2010), this object exhibited a redder- 
when-brightcr trend when it was fainter than V ~ 15, 
while a bluer-when-brighter trend was observed when 
it was brighter. This feature can be confirmed in fig- 
ure 7. A redder-when-brighter trend was detected in the 
whole data because the object had almost stayed in the 
faint state during our observation period. As well as 
3C 454.3, PKS 1510-089 and PG 1553+113 also exhib- 
ited redder-when-brighter trends in their faint states (see, 
figures 21 and 28). No significant correlation of the flux 
and color was detected in those two objects when their 
whole data were used. Bluer-when-brighter trends were 
significantly detected in 3C 454.3 and PKS 1510-089 in 
their bright states (V < 15.0 for 3C 454.3 and < 15.5 for 
PKS 1510-089). 

It is widely accepted that the rcddcr-whcn-brightcr 
trend in the faint state is attributed to a strong contri- 
bution of thermal emission from an accretion disk (also 
see, Smith et al. 1986; Hagcn-Thorn, Yakovleva 1994, 
Pian et al. 1999). The contribution of the thermal disk 
emission would definitely be small when the flux from the 
jet increases. Hence, the bluer-when-brighter trends in 
the bright states of 3C 454.3 and PKS 1510-089 prob- 
ably have the same origin as those in ordinary objects. 
Including 3C 454.3 and PKS 1510-089, the number of 
objects that showed a blucr-when-brighter trend is 25. 

PG 1553+113 exhibited a significant redder-when- 
brighter trend in its faint state (V > 13.8), as well 



10 



Y. Ikejiri et al. 



[Vol. 



13.8 

14 

14.2 

14.4 

14.6 
1.5 
1.6 
1.7 
1.8 
1.9 
2 

13.4 
13.6 
13.8 
14 
14.2 
14.4 
14.6 



I I I 

• 

• * 




— I 1 h 

: ... I. 


1 h- 














'A 





4700 



4800 4900 5000 
JD-2450000 



5100 



~i 1 1 1 r 



_I I I I L_ 



1.4 1.5 1.6 



1.7 

V-J 



1.8 1.9 




13.5 
14 

14.5 
15 

15.5 
16 

16.5 



4600 4700 4800 4900 5000 5100 5200 
JD-2450000 



.•'.••.7«'r:." 

_l I I I I L_ 



1.2 1.4 1.6 



1.8 2 
V-J 



2.2 2.4 2.6 



4650 4700 4750 4800 4850 4900 4950 
JD-2450000 





16 
16.5 

17 

17.5 

18 
2 

2.2 

2.4 

2.6 



1 

* 


^ i i 

h 

8 

CD 

1 - 

— f h— 




i i 



4700 
14.5 

15 
15.5 

16 

16.5 
1.6 



4800 4900 5000 
JD-2450000 



5100 



4600 4700 4800 4900 5000 5100 5200 
JD-2450000 



4600 



4800 5000 
JD-2450000 



~i 1 1 r 



.v. 



_] I I L_ 



5200 



1.7 



1.8 1.9 
V-J 



2.1 




Fig. 7. Examples of the V-band light curve (top), V — J color variation (middle), and their color-magnitude diagram (bottom) of 
blazars. Panel a: 3C 371, b: 3C 454.3, c: PG 1553+113, d: PKS 0048-097, e: H 1722+119, and f: PKS 1502+106. In the lower 
three panels, different symbols represent observations taken in different periods of time. 
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as 3C 454.3 and PKS 1510-089. No bluer-when- 
brighter trend was significantly detected in its bright 
state. Panel (c) of figure 7 shows the color-magnitude di- 
agram of this object. The redder-when-brightcr trend in 
the faint state, however, possibly weakens above V = 13.8. 
This weakening of the redder-when-brighter trend may 
suggest a transition of a bluer-when-brighter trend in a 
bright state of PG 1553+113. 

There are seven blazars that showed no significant cor- 
relation between the V-band magnitude and the V — J 
color, in addition to PKS 1510-089 and PG 1553+113. 
The numbers of observations are quite small for two ob- 
jects: PKS 0215+015 and S4 0954+65. Only 6 and 5 
observations were made, respectively. They are too small 
to conclude anything about the correlation. There are 
five blazars that was observed for > 10 d and showed 
no significant correlation of the light curve and color: 
QSO 0454-234, MisV1436, S5 1803+784, PKS 0048-097, 
and H 1722+119. 

Even in these five objects, bluer-when-brighter trends 
were occasionally seen in a short period of time. Panel (d) 
of figure 7 shows the variation in the magnitude and 
the color for PKS 0048-097. Our observation of 
PKS 0048-097 covered two seasons in 2008 (JD 2454741- 
2454851) and 2009 (JD 2455068-2455097). A significant 
bluer-when-brightcr trend was detected in the 2008 ob- 
servation shown in the filled circles in panel (d). A pos- 
sible bluer-when-brighter trend can be seen also in the 
2009 observation shown in the open circles, while the 
number of observations is too small to conclude it. It is 
important to note that, in the color-magnitude diagram, 
the 2009 data apparently shows a different bluer-when- 
brighter slope from the 2008 data. This is the reason why 
no significant correlation was detected when we performed 
a correlation analysis using all data. 

A similar behavior can also be seen in H 1722+119 
[panel e of figure 7]. We divided our observa- 
tion period into four periods of time (JD 2454659- 
2454720, JD 2454747-2454760, JD 2454908-2455001, and 
JD 2455079-2455121), and show them by different sym- 
bols in panel (e). Bluer-when-brighter trends are appar- 
ently seen in the first (filled circles), second (open cir- 
cles), and third (filled squares) epochs, while the number 
of observations is small. A sign of bluer-when-brighter 
trends can also be seen in QSO 0454-234, MisV 1436, and 
S5 1803+784 in their bright states (see figure 22, 25, and 
27, respectively). QSO 0454—234 exhibited a significant 
bluer-when-brightcr trend in its bright state (V < 16.0), 
while the trend is not statistically significant in MisV 1436 
and S5 1803+784. 

We found that PKS 1502+106 also exhibited a flux- 
color behavior similar to those of PKS 0048—097 and 
H 1722+119, although a redder-when-brighter trend was 
detected in this object in its whole data from 2008 to 2010. 
Panel (f) of figure 7 shows the variations in the magnitude 
and color of PKS 1502+106. We calculated correlation co- 
efficients for each year, and found significant bluer-when- 
brightcr trends in both the 2008 and 2009 data sets. In 
2010, this object remained faint (V > 17.2) and showed no 



80 

70 

60 

50 

40 

30 

20 

10 


5 

-5 



1 1 


i i 

.-+" 

+ + "" 


- 


+ &c 




+ W 
4+ .-••+ + 


+ y 
+++' 


+ 


++++.+' 




H 1 — 

" + +* ++ + \ + 


l-^f 1 

t+ 


+ 

i i 


++ + + + + + + + " 
+ + + 
, + 



10 20 30 40 

V-band flux (x10" 12 erg s" 1 cm" 2 ) 



50 



Fig. 8. Upper panel: J- and V-band flux variations observed 
in PKS 1749+096. The lower panel: Residuals of the J-band 
flux from the best-fitted linear model. 

significant correlation between the light curve and color. 

The sample reduces to 32 blazars once we exclude 10 
objects, in which the number of V — J observations is 
< 10. In these 10 objects, observations are too few to con- 
clude whether there is a significant correlation between the 
light curve and color. Among those 32 blazars, in sum- 
mary, 23 blazars showed significant bluer-when-brightcr 
trends in their total data. Five blazars showed temporary 
significant bluer-when-brightcr trends in short periods of 
time (QSO 0454-234, PKS 1510-089, PKS 1502+106, 
3C 454.3, and PKS 0048-097). Then, we conclude that 
a total of 28 blazars, corresponding to 88 % of the well- 
observed sample, had the bluer-when-brighter trend. Four 
blazars (PG 1553+113, H 1722+119, S5 1803+784, and 
MisV 1436) possibly showed a sign of bluer- when-brighter 
trends, although the trends were not statistically signifi- 
cant in our data. 

3.2.3. J- and V-band Flux Diagram 

Most blazars exhibited the bluer-when-brighter trend, 
as mentioned in the last subsection. Two scenarios can 
be considered to generate the bluer-when-brighter trend. 
The first scenario is that the variation component, itself, 
had a bluer-when-brighter trend. The second scenario is 
that the color of the variation component was constant 
with time, while an underlying component was redder 
than the flare component. We can evaluate those two 
scenarios using the Fj-Fy diagrams, where Fj and Fy 
denote the J- and V-band fluxes, respectively. Examples 
of the Fj-Fy diagrams are presented in figures 8, 9, and 
10. In these figures, achromatic variations are described 
in a form of Fj = cFy, where c is a constant. 

Fj can be described with a linear function of Fy in most 
objects. A typical example is PKS 1749+096, as shown in 
figure 8. The linear relationship holds over one order of 
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Fig. 9. As for figure 8, but for PKS 1510-089. 
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Fig. 10. Top panel: As for figure 8, but for 3C 454.3. The 
filled and open circles represent the observations in 2008 and 
2009, respectively. The middle and bottom panels: Residuals 
of the 2008 and 2009 data, respectively. 



magnitude of the flux. We confirmed that linear regression 
models for Fj and Fy yielded significant positive intersec- 
tions of Fj at Fy — in all bluer-when-brighter objects. 
These results support that there was an underlying com- 
ponent which was redder than the variation components 
and the color of the variation components was roughly 
constant. 

We note that there are two possible exceptions of our 
28 blucr-when-brighter blazars in terms of the Fj-Fy re- 
lation: PKS 1510-089 and 3C 454.3. Their Fj-F v di- 
agrams are shown in figures 9 and 10, respectively. In 
both objects, we can see a sign of systematic deviations 
from their best-fitted linear models; the slope is steep in a 
low flux regime, but shallow in a high flux one. The 2008 
and 2009 observations of 3C 454.3 are represented by the 
filled and open circles, respectively. PKS 1510—089 and 
3C 454.3 had a common feature that they showed redder- 
when-brightcr trends when they were faint, as mentioned 
in sub-subsection 3.2.2. The low-flux regimes in the 
Fj-Fy diagram correspond to the reddcr-when-brighter 
phase. 

Next we focus on the behavior of short-term flares on 
time scales of days-weeks in the Fj-Fy diagram. In sub- 
subsection 3.2.1, we describe our search for spectral hys- 
teresis associated with flares through time-lags between 
different band fluxes, or between the flux and the color. 
This approach could be effective if the observed colors 
directly represent the color of the variation components. 
However, as mentioned above, the observed color is prob- 
ably a composition of an underlying red component and 
variable ones. It could be better to search for a sign of 
hysteresis using the Fj-Fy diagram. 

Figure 1 1 shows the Fj-Fy diagram (lower panels) and 
the V-band light curves (upper panels) for three short- 
term flares (flare a, b, and c) observed in PKS 1749+096. 
Flares a, b, and c are shown in left, middle, and right 
panels, respectively. In flare a, the rising phase appears 
to draw a different path from the decay phase in the Fj- 
Fy diagram. As a result, we can see hysteresis in this 
diagram. The object appears to draw a "loop" in the Fj- 
Fy diagram, moving in the anti-clockwise direction in the 
loop during the flare. This behavior means that the rising 
phase of the flare was bluer than the decay phase around 
the flare maximum. We can also see similar hystereses in 
flares b and c, while the hysteresis in flare c may be just an 
erroneous one because a detailed behavior was overlooked 
around the peak of flare c. 

We searched for such hystereses for all short-term flares 
that wc observed. It is expected that a hysteresis of a 
small flare in the Fj-Fy diagram is disturbed by the su- 
perposition of other small flares. Hence, we again used 
the 17 prominent flare events that are selected in sub- 
subsection 3.2.1. Then, we calculated the slopes in the 
Fj-Fy diagram for rise and decay phases of those flares. 
The slopes were defined as: 

(AFj/AF v ) Tise = (Fj tP - F JiP _i)/(Fy, p - F v>P -i) (4) 

(AFj/ AFy)dccay = (Fj, p - Fj,p+i) I {Fv,p — -FV,p+l),(5) 

where Fjry) tV and Fj(y^ p ±i denote the J-(^-)band fluxes 
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Fig. 11. V-band light curves (upper panels) and the Fj-Fy diagrams (lower panels) of PKS 1749+096. The rise and decay phases 
of three flares (flare a, b, and c) are represented by the solid and dashed lines in each panel, respectively. 




AFj/AF v (Rise) 



Fig. 12. Slope in the Fj-Fy diagram of the rise and de- 
cay phases around the maxima of short- term flares. The 
filled circles and crosses represent the flares which were 
selected with peak-to-peak amplitudes of > 1.0 mag and 
> 0.1 mag, respectively. The dashed line represents 
(AFj/AF v ) lise = (AFj/AF v ) 

decay • 

at the flare maximum and its neighboring observations. 
The case of (AFj / AF v ) Tisc = (AFj / AFy) deca , y means 
that there is no hysteresis. The case of (AFj / AFy) lisc > 
(AFj/ AFy) decay indicates that a rising phase is bluer 
than a decay phase, as can be seen in figure 11. The 
case of (Ai 7 'j/AFy) r i S0 < (AFj / AFy) decay indicates that 
a rising phase is redder than a decay phase. 

Figure 12 shows AFj/ AFy of the rise and decay phases 
of the 17 flares, as indicated by the filled circles. There are 
11 flares showing (AFj / AF v ) Hsc > (AFj / AF v ) deca , y , five 
flares showing (AFj / AF v ) risc < (AFj / AF v ) dccay , and 



one flare showing no significant hysteresis. This result 
suggests that the spectral hysteresis of (AFj / AFy) Tisc > 
(AFj I AFy) decay was more preferable for the prominent 
short-term flares. The crosses in the figure represent 
AFj I AFy for small flares selected with a peak-to-peak 
amplitude of 0.1-1.0. The fraction of these small flares 
showing (AFj / AFy) risc > (AFj / AF v ) dccay was compa- 
rable to that showing (AFj/A_Fv) r i so < (AFj/AFy^ccay 
This is probably due to the superposition of multiple flares 
that disturb the hysteresis pattern of each flare. 

The spectral hysteresis was detected in the sample of 
prominent flares whose peaks were recorded within 5 d 
after and before the neighboring observations (see sub- 
subsection 3.2.1). Together with the time-lag analysis de- 
scribed in subsection 3.2.1, the time-scale of the hysteresis 
events that we detected is, hence, less than a few days. 
This is comparable with or possibly shorter-than-average 
sampling rate of our monitoring. It is possible that short- 
term hysteresis events could be overlooked in the flares of 

(AFj/AFy) riso < (AFj / AFy) 

decay ■ 

3.3. Light Curve and Polarization 

3.3.1. Time lag between variations in the light curve and 
the polarization degree 
We searched for time lags between the light curve and 
the polarization degree (PD) using the DCF method, as 
described in sub-subsection 3.2.1. Figure 13 presents an 
example of DCF, which was calculated with all data of 
the V-band flux and PD of BL Lac. As shown in this fig- 
ure, we found no significant time-lag between them. We 
calculated DCF for all blazars that we observed, and con- 
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Fig. 13. DCF calculated with the V-band light curve and 
the polarization degree of BL Lac. 

firmed that no objects exhibited a significant time-lag be- 
tween the light curve and PD. As mentioned in the next 
subsection, the correlation, itself, was weak between the 
flux and PD. 

We selected 17 prominent flare events in sub- 
subsection 3.2.1. The peak timings of their light curves 
were coincident with those of PD, or there was no corre- 
lation between them at all, except for the flares observed 
in PKS 1510-089, AO 0235+164, and PKS 1502+106. 
Sasada et al. (2011) report that the increases in PD were 
associated with prominent flares in PKS 1510—089 and 
AO 0235+164 using the same data as that in the present 
paper. The PD peaks preceded the flux peaks by 1 d in 
PKS 1510-089 and 4 d in AO 0235+164. On the other 
hand, figure 14 shows the cases that the PD maximum 
lagged the flux one. The upper panel shows the light curve 
and PD in PKS 1502+106. A large flare having double 
peaks can be seen in JD 2454940-2454970. The double- 
peaked feature was also seen in the variation in PD, while 
both peaks of PD lagged behind those of the light curve 
by 3 d. Similar events were also found in 3C 454.3. A 
large flare in ~JD 2455170 was followed by a subsequent 
increase in PD, as shown in the lower panel of figure 14. 
The lag was about 5 d. All of those time-lags were not 
universal in each object; there were flares showing no or 
shorter time-lags in all objects. 

Thus, there was no universal time-lag between the flux 
and PD. As well as the flux-color correlation, the flux- 
PD correlation is also discussed without respect to time 
lags in the next subsection. 

3.3.2. Correlation between the V-band light curve and 
the polarization degree 
We report on the correlation of the flux and PD in this 
subsection. We present examples of the light curves and 
variations in PD in figure 15. Panel (a) of figure 15 shows 
those of AO 0235+164. PD of this object became higher 
when it was brighter. As described in subsection 3.2.2, we 
calculated the Pearson correlation coefficients with the V- 
band flux and PD, and tested them using the Student's 
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Fig. 14. V-band light curves and PD variations of flares in 
PKS 1502+106 (upper panel) and 3C 454.3 (lower panel). 
The filled and open circles represent the light curve and PD, 
respectively. 

i-test. The correlation coefficients are presented as r po i. in 
table 2. We identified that 10 blazars exhibited significant 
positive correlations between the V-band flux and PD. 

The flux-PZ) correlation was not simple in several ob- 
jects, even for r po i. > 0. For example, panel (b) of figure 15 
shows the V-band light curve and PD of MisV 1436. 
This is one of objects that showed a significant posi- 
tive correlation between the flux and PD. However, 
PD in JD 2454817-2454901 were much lower than those 
in JD 2455034-2455085, although the V-band magni- 
tudes in these periods were comparable. Similar behavior 
was also observed in PKS 1749+096. Significant posi- 
tive correlations between the flux and PD were detected 
when we analyzed the flares between JD 2454720-2454760 
and JD 2454894-2455084 in PKS 1749+096 (figure 23). 
However, no significant correlation was detected when we 
used all data. This is due to the presence of flares with- 
out an increase in PD. Hence, the correlation was not 
universal in those objects. 

We identified that four blazars (OJ 49, ON 231, 3C 66A, 
and ON 325) showed a significant negative correlation of 
the flux and PD, namely r po i, < 0. Panel (c) of figure 15 
shows the light curve and PD of OJ 49. PD were lower 
when the object was brighter. In addition to those four 
blazars, OJ 287 and BL Lac also exhibited flares in which 
negative correlations between the flux and PD were ob- 
served: A flare of OJ 287 in JD 2454612-2454983 and a 
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Fig. 15. Examples of the time variation in the V-band magnitude (top) and polarization degree (middle), and the flux-polarization 
degree diagram (bottom). Panel a: AO 0235+164, b: MisV1436, and c: OJ 49. 



flare of BL Lac in JD 2454612-2454864 (figures 23 and 
26). 

The sample reduces to 33 blazars once we exclude the 
objects in which available polarization data were less than 
10. The numbers of objects with significant positive and 
negative correlation correspond to 30 % and 12 % of the 
sample, respectively. 

In addition, we calculated the correlation coefficients 
between the V — J color and PD. The results arc shown 
in table 2 as r co i._ po i.. A positive (or negative) r co i._ po i. 
means that an object becomes redder (or bluer) when PD 
increases (or decreases). The numbers of objects with 
significant positive and negative correlation, and non- 
significant correlation are, 2, 3, and 29, respectively. The 
correlation between V — J and PD was quite weak. 
3.3.3. Rotation of polarization angle 

Possible rotations of polarization angle (PA) have been 
reported and discussed both in optical and radio observa- 
tions. However, a part of them was so sparsely sampled 
that it can be explained by a result of a random walk in 
the Stokes QU plane (e.g. Aller et al. 1981; Jones et al. 
1985; Kikuchi et al. 1988). Hence, polarimetric observa- 
tions with high observation density are required to estab- 
lish real rotation events of polarization. Recently, smooth 
rotations of optical polarization have been reported in sev- 
eral blazars based on long-term polarimetric observations 
(Marscher et al. 2008; Sasada et al. 2010; Marscher et al. 
2010; Abdo et al. 2010c; Jorstad et al. 2010). Our po- 
larimetric observation allowed us to perform a systematic 
search for such rotations. 

A search for polarization rotation events is, however, 
not straightforward. Two issues should be taken into ac- 
count; one is an ambiguity of 180° in PA, and the other is 
the presence of multiple polarization components. First, 
we resolved the 180° ambiguity by assuming that objects 



favor gradual changes in PA rather than rapid changes. 
In other words, a difference in PA between neighboring 
two observations should be smaller than 90°. We added 
-180° (or +180°) to the observed PA in the case that 
PA n > PA n -! +90° ( or PA„ < PA n _ x -90°), where 
PAi is the «-th observation. The corrected PA by this 
method is shown in the left-bottom panel in the figures in 
Appendix. 

Second, there are blazars in which a stationary or long- 
term variation component of polarization has been pro- 
posed (e.g. Moore et al. 1982; Villforth et al. 2010). 
Polarization rotations in such objects should be discussed 
for each component, i.e., long- or short-term ones. In this 
paper, we define the average QU as a stationary polariza- 
tion component. The averages of Q and U were calculated 
iteratively; we first calculated an average of Q (or U) us- 
ing all data, and after discarding outliers (> 3<r) from the 
sample, re-calculated it. The iteration continued until no 
outliers remained. 

Even after those two-types of corrections were per- 
formed, we could not find any clear universal correlations 
between the light curve and any sets of corrected PA. 
However, we found several noteworthy objects or flares 
in which possible rotations of polarization were seen. We 
describe them below. 

Figure 16 shows the temporal variation in the ob- 
served and corrected PA and the QU distribution of 
PKS 1510-089. The observed PA, shown in panel (b), ap- 
parently varies randomly. However, as shown in panel (c), 
we can see an increasing trend in the 180°-ambiguity cor- 
rected PA during an active phase of the object. This rota- 
tion event has also been reported in Marscher et al. (2010), 
and confirmed in Sasada et al. (2011). The rotation event 
continued for ~ 90 d and its amplitude in PA was over 
700°. The long duration and large amplitude of the ro- 
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Fig. 16. Variation in PA of PKS 1510-089 observed in 2009. 
Panels (a) the V-band light curve, (b) observed PA, (c) cor- 
rected PA for the 180° ambiguity, (d) PA measured from the 
average (Q,U), and (e) the QU plane around the origin of 
QU. In panel (e), the average (Q,U) is indicated by the filled 
circle. 



tation event are definitely difficult to be explained by the 
result of a random walk in the QU plane (e.g. Moore et al. 
1982). 

Similar rotation events having a long duration and a 
large amplitude were seen only in two objects: 3C 454.3 
(figure 25) and PKS 1749+096 (figure 23). All of those 
rotation events have a common feature, that they were 
associated with active phases or flares. This implies that 
the rotation events were related to the variation in the 
magnetic field direction in the emitting region of the 
flares. Another possible rotation event can be seen in 
S5 0716+714 (figure 29). The change in PA of this ob- 
ject was, however, not as smooth as observed in the above 
three objects, but rather had a step-like profile. It contin- 
ued for a very long time of ~ 400 d. 

It should be noted that the 180°-ambiguity correc- 
tion is sensitive to the sampling density of observations. 
Furthermore, the results should be considered carefully 
if there are multiple polarization components. The be- 
havior of PA in PKS 1510—089 apparently changes if 
the origin of (Q,U) is shifted to their observed aver- 
ages, as shown in panel (d) of figure 16. The increas- 
ing trend in PA is terminated earlier than that shown in 
panel (c). The average (Q,U) was (3.9 x 10" 14 ,2.6x 10~ 14 ) 
(in erg s _1 cm -1 Hz -1 ), indicated by the filled circle in the 
QU plane (panel e). This average (Q,U) corresponds to a 
PD of ~ 0.6% at V = 16.0. This result demonstrates that 
the 180°-ambiguity correction could be affected by an- 
other polarization component, even if it has a quite small 
contribution to the total flux. The polarization rotation in 
PKS 1510—089 probably started at the onset of the active 
phase. It was then interrupted by a sharp spike-like flare 
in ^JD 2454960. After this spike, it is unclear whether 
the rotation continued or not. 

Another noteworthy case is 3C 66 A. Figure 17 is 
the same as figure 16, but for 3C 66A. Our 180°- 
ambiguity correction yields a short-term rotation event 
between JD 2454840-2454850, as can be seen in panel (c). 
However, this feature cannot be seen in panel (d), in which 
PA are measured from the observed average of (Q,U). As 
shown in the QU plane in panel (e), the observed QU are 
distributed not around the origin of (Q,U), but around 
the average point indicated by the filled circle. Hence, 
this observation suggests that the observed polarization 
in 3C 66A can be interpreted as a composition of two 
polarization components, namely, long- and short-term 
variation components. If this is the case, the apparent 
PA variation in panel (c) means no real variation in the 
magnetic field direction in the emitting area. A long-term 
component should be subtracted from observed polariza- 
tion when the rotation of polarization is discussed for ob- 
jects like 3C 66A. 

Similar cases to 3C 66A can be found in 
RX J1542.8+612 (figure 27) and QSO 0454-234 
(figure 24); short-term rotation events appeared in 
their corrected PA. Furthermore, their average (Q, U) 
were significantly deviated from their origin of (Q,U). 
Marschcr et al. (2008) have reported the detection of a 
polarization rotation in BL Lac. The rotation event was 
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Fig. 17. Variation in PA of 3C 66A. Panels (a) the V-band 
light curve, (b) observed PA, (c) corrected PA for the 180° 
ambiguity, (d) PA measured from the average (Q,U), and 
(e) the QU plane in which all observations are included. In 
panel (e), the average (Q,U) is indicated by the filled circle. 



well sampled, continuing for about one week and having 
an amplitude in PA of 240°. This object has, however, 
been proposed to have two polarization components. It is 
known that the observed average of (Q,U) is significantly 
deviate from the origin of (Q, U) (Hagen-Thorn et al. 
2008). Our observation confirmed it, as can be seen 
in figure 26. If short-term variation components are 
superimposed on the stationary or long-term polarization 
component, the model for the rotation event reported in 
Marscher et al. (2008) should be reconsidered. 

Finally, we show the case for PKS 1502+106 in fig- 
ure 18. A rotation of polarization is not seen in the 
observed and 180°-ambiguity corrected PA [panels (b) 
and (c)], but can be seen in the origin-shifted PA [panel 
(d)]. The validity of this correction is, however, brought 
into question. The observed (Q, U) is expected to be dis- 
tributed centered around the stationary component in the 
two-component picture. However, there is only few points 
at the vicinity of the average (Q,U), as shown in panel (e). 
This result is inconsistent with the picture of the two- 
component model. This example tells us that the shift of 
the origin of (Q,U) could generate false variations in PA 
in some cases. 

In summarizing this subsection, we offer several cau- 
tions to avoid misinterpretation of observed PA variations 
in terms of the polarization rotations. First, it is required 
that a rotation episode is well sampled when it is proposed 
to be a real one. A rotation episode indicated by only a 
few data points, cannot be distinguished from results of a 
random walk in the QU plane. Second, a rotation event 
would more likely to be real when it has a larger vari- 
ation amplitude in PA. An apparent rotation having a 
small amplitude of < 180° can be readily produced by a 
composition of a few flares having different PA. Third, 
the appearance of a rotation event should not be changed 
by a small shift of the origin of (Q,U). As mentioned 
above, a part of the rotation event in PKS 1510—089 is 
quite sensitive to a small change of the origin of (Q,U). 
Finally, the temporal behavior in the QU plane should be 
carefully checked not only during the rotation event, but 
also before and after it. If the observed average of (Q,U) 
is significantly deviated from the origin, a stationary or 
long-term variation component in polarization could be 
present. Then, polarization variations should be discussed 
after the long-term component is subtracted from the ob- 
servation. 

3.4- The Vpeak and luminosity dependence of variability 

In this subsection, we report on the dependence of the 
degree of variability in the flux, color, and PD on t'peak- 
The ^poak is related to the luminosity of the blazars in our 
sample (see subsection 3.1). We note that variation am- 
plitudes may be underestimated because blazars generally 
have a variation time-scale longer than our observation 
period. To minimize this observation bias, we calculated 
the slope in the color-magnitude diagram, dVj d(V — J), as 
well as ordinary peak-to-peak amplitudes. We estimated 
the slope, dV/d(V — J), by fitting a linear function to the 
observed V and V — J. dV/d(V — J) can be considered as 
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Fig. 19. The fp Ca ]j and My distribution of the degree of variability. From top to bottom, the panels show dV/d(V — J), the 
peak-to-peak variation amplitudes of the V-band flux (AV), the V - J color [A(V- J)], and the \/-band PD (APD). The filled 
circles, triangles, and open circles represent BL Lac objects, FSRQs, and poorly observed objects in which the number of observations 
are less than 20, respectively. 
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Fig. 18. As for figure 17, but for PKS 1502+106. 



an index of the degree of the flux variability against the 
color variation. 

Figure 19 shows the dependence of the degree of varia- 
tions on i-'poak and the average absolute magnitude, My, 
in the left and right panels, respectively. From top to bot- 
tom, the panels show dV/ d(V — J), the peak-to-peak vari- 
ation amplitudes of the V-band flux (AV), the V — J color 
[A(V — J)], and the V-band PD (APD). In the case of 
3C 454.3 and PKS 1510-089, dV/d(V — J) were calculated 
using the data in their bright states (V < 15 for 3C 454.3 
and < 15.5 for PKS 1510-089). This is for eliminating the 
redder-when-brighter stage (see, sub-subsection 3.2.2). 

All panels in figure 19 suggest that the degrees of vari- 
ation were low in high ^ poa k, or faint objects. They can 
be high in low ^ pca k, or bright objects. The poorly ob- 
served objects represented by the open circles tend to be 
distributed in low variability regions compared with the 
others. This is naturally interpreted as undercstimations 
of the variation amplitudes in those objects due to the 
short period of observations. Several FSRQs, particularly, 
3C 273 and QSO 0454—234, exhibit a small variation am- 
plitude, although they have a quite low v P eak- This atyp- 
ical feature can be understood by large contamination of 
the thermal emission from accretion disks. Those invari- 
able component probably made the variability of the jet 
component small, apparently. Thus, the variation ampli- 
tudes of poorly observed objects (the open circles) and 
FSRQs (the filled triangles) can be considered as lower 
limits of the amplitudes. Then, we can see good corre- 
lations of dV/d{V - J), AV, A(V - J) , and APD with 
fpeak and My. Those characteristics of the degree of vari- 
ation suggest that the variability is stronger in a low i/ pea k 
object, and hence in the emission from higher energy elec- 
trons. 

The fact that the variability in the flux and color is 
larger in lower ^ poa k objects has been reported in previ- 
ous studies (e.g. Ulrich et al. 1997). The present study 
confirmed this feature with a large sample. In addition, 
the same feature can be also seen in the variability of PD. 
This implies that the flares in blazars are generally asso- 
ciated with the variation of PD (see also, Sasada et al. 
2011). 

4. Discussion 

4-1. "Bluer- When- Brighter", as a Universal Aspect of 
Blazars 

Our V- and J-band photometric observation showed 
that the "bluer-when-brighter" trend was detected in a 
high fraction of blazars we observed: 28/32 well-observed 
sources. The fraction of bluer-when-brighter objects is 
quite high compared with those reported in previous stud- 
ies: Gu et al. (2006) investigated variations in V — I of 
three FSRQs and five BL Lac objects. They found that 
three BL Lac objects exhibited significant bluer-when- 
brighter trends, while two FSRQs exhibited redder-when- 
brighter trends. Rani et al. (2010) performed an anal- 
ogous analysis for six FSRQs and six BL Lac objects. 
They found that three BL Lac objects exhibited signif- 
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icant bluer-when-brighter trends, and four FSRQs exhib- 
ited rcddcr-when-brighter trends. Those studies thus pro- 
vide lower fractions of bluer-when-brighter objects than 
our results. On the other hand, several studies have in- 
dicated that the bluer-when-brighter trend is common in 
BL Lac objects: D'Amicis et al. (2002) have investigated 
variations in the optical spectral index of eight BL Lac 
objects, and found that all objects showed bluer- when- 
brighter trends. Fiorucci et al. (2004) report on their 
photometric study of 37 blazars, and propose that bluer- 
whcn-brightcr trends are more prominent in objects show- 
ing stronger variability. In those two studies, no evalua- 
tion of the significance of the correlation has been given. 
Vagnetti et al. (2003) also analyzed the spectral slope vari- 
ability for eight BL Lac objects, and reported that all 
objects showed bluer- when-brighter trends. However, the 
sampling rate of their data was ~ 25 d, and it is too large 
to study variability on a time scale of days. Thus, our 
present study, for the first time, yields an unambiguous 
confirmation that the bluer-when-brighter trend is a com- 
mon feature in variations on time scales of days-months 
in blazars. 

We propose that there are two factors that can disturb 
the general bluer- when-brighter trend. First, a strong 
contribution of a blue thermal emission from the accre- 
tion disk leads to a redder-when-brighter trend. A typical 
case is 3C 454.3, as shown in panel (b) of figure 7. The 
observed color variation suggests that the redder-when- 
brightcr trend was observed only when the objects were 
faint. This color behavior is consistent with the scenario 
that the contribution of the thermal emission from the ac- 
cretion disk becomes strong when the synchrotron emis- 
sion from the jet weakened (Villata et al. 2006a). Even 
in the redder-when-brightcr objects, our observation re- 
vealed that the bluer-when-brighter trend appeared when 
they were bright. 

Second, the bluer-when-brighter trend could be hidden 
if there are multiple variation components having differ- 
ent bluer-when-brighter sequences. We see the color be- 
havior of long-term variation components, taking the case 
of PKS 1502+106 as an example. Panel (f) of figure 7 
shows the light curve, color, and color-magnitude diagram 
of PKS 1502+106. A bluer-when-brighter trend was not 
detected in the whole data of this object, but significantly 
detected in each year, as reported in sub-subsection 3.2.2. 
Figure 20 shows the Fj~F v diagram of PKS 1502+106. 
The filled circles and crosses represent data obtained in 
2008 and 2009, respectively. As can be seen in this figure, 
Fj can be described with a linear function of Fy both 
in the 2008 and 2009 data. The dashed lines in the fig- 
ure represents the best-fitted linear regression models for 
each year. The slope of the 2008 data in the Fj-Fy dia- 
gram is clearly different from that of 2009; a steep slope 
in the 2009 data compared with the 2008 one. As can be 
seen from figure 7, short-term flares are apparently super- 
imposed on a long-term outburst component in the 2009 
light curve. The long-term component reached the maxi- 
mum at the beginning of our 2009 observation, and then 
started a gradual decay through that year. Hence, the 
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Fig. 20. Upper panel: J- and V-band flux observed in 
PKS 1502+106. The filled circles and crosses represent the 
observations in 2008 and 2009, respectively. The dashed lines 
indicate the best-fitted linear regression models for each year. 
The thick solid lines represent rising phases of large flares. 
Lower panel: Residuals of the J-band flux from the best-fit- 
ted model for the 2009 data. The dashed lines represent the 
slope of the best-fitted model for the 2008 observations. 

general trend seen in the Fj-Fy diagram of the 2009 data 
actually originates from the long-term outburst compo- 
nent. In contrast, the object almost remained fainter than 
V = 17 in 2008, while several short-term flares brighter 
than V = 17 were observed. There was no prominent out- 
burst component in 2008. 

It is interesting to note that, in the 2009 season, the Fj— 
Fy slopes of the early phase of the short-term flares were 
possibly analogous to those in 2008. The thick solid lines 
in figure 20 represent the rising phases of the prominent 
short-term flares peaked at JD 2454951 and 2454960. The 
lower panel shows the residuals of Fj from the best-fitted 
model for the 2009 data. The dashed line in the lower 
panel indicates the Fj-Fy slope of the 2008 data. The 
figure suggests that there were stages of shallow Fj-Fy 
slopes during the early rising phases of the flares in 2009. 
After these stages, the slope appears to become very steep 
near the flare peaks. This may be a sign of spectral hys- 
teresis, which was frequently observed in short-term flares, 
as mentioned in sub-subsection 3.2.3. The early shallow 



No. ] 

Fj-Fy slopes look close to the slope of the short-term 
flare in 2008. 

Those observations give us an idea that the color behav- 
ior of the short-term flares in 2008 was the same as that in 
2009. In addition, the color of the outburst component in 
2009 was probably redder than the short-term flares. As 
well as PKS 1502+106, PKS 0048-097 and H 1722+119 
clearly showed different color behavior in each year, as 
shown in the lower panels of figure 7. Those phenomena 
may also be due to the presence of a long-term compo- 
nent, whose color was different from that of short-term 
flares, and moreover, variable with time. This idea is sup- 
ported by the fact that correlation coefficients were low 
even in the bluer- when-brighter sources; as shown in ta- 
ble 2, r co \ was less than 0.70 in more than half of objects. 
The low correlation-coefficients can be interpreted by the 
scenario that the observed color is a composition of short- 
and long-term components that have different bluer-when- 
brighter sequences. The two-component scenario has also 
been proposed to explain the color behavior of BL Lac 
(Villata et al. 2002; Villata ct al. 2004). 

The above two factors, namely the presence of the 
thermal disk emission and multiple synchrotron varia- 
tion components, can disturb the general bluer-when- 
brighter trend in blazar variability. Our observation over- 
comes those factors by covering a wide range of wave- 
lengths, namely, the simultaneous optical and NIR ob- 
servations and a high sampling rate. Those two advan- 
tages of our observation led to the high fraction of bluer- 
when-brighter objects compared with previous studies. In 
other words, one should carefully discuss apparent redder- 
when-brighter trends or achromatic variations observed in 
blazars, while taking into account the above two factors 
that can disturb the universal bluer-when-brightcr trend. 

It has been suspected that the bluer-when-brighter 
trend is common in BL Lac objects, but not in FSRQs 
because of a large contribution of thermal emission in 
FSRQs (Gu et al. 2006; Rani et al. 2010). Our objects that 
were observed more than 10 times include seven FSRQs 
and seven LBLs. The correlation analysis with all data 
sets showed that significant bluer-when-brightcr trends 
were detected in two FSRQs and six LBLs (see table 2). 
This result confirms the previous studies, which suggest 
that the color-flux correlation of FSRQs is different from 
that of BL Lac objects. However, the bluer-when-brighter 
trend was significantly detected even in the other four 
FSRQs when they were bright (3C 454.3, PKS 1510-089, 
PKS 1502+106, and QSO 0454-234), as mentioned in 
sub-subsection 3.2.2. Hence, the jet emission in FSRQs 
presumably has the same color behavior as that in BL Lac 
objects. 

The relationships of Fj and F v of 3C 454.3 and 
PKS 1510—089 suggest that their variation components 
did not have a constant color, but a blueing trend with an 
increase in flux (see figure 9 and 10). Alternatively, it can 
be interpreted with the scenario that two variation com- 
ponents had different time-scales and colors, as proposed 
above. A red and long-term outburst component would 
be dominant in their low flux regime. The redder-when- 
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brighter trend was observed probably because an under- 
lying source was the blue thermal emission. The emission 
from short-term flares would be dominant in the high- 
flux regime, resulting in the standard bluer-when-brightcr 
trend. This scenario is supported by the light curve struc- 
ture; short-term flares, making the objects brightest, were 
superimposed on long-term outburst components (see fig- 
ures 25 and 23). 

In general, variations on time-scales of months-years 
are expected to have larger amplitudes than those of days- 
weeks in blazars (Hufnagcl, Brcgman 1992; Takahashi 
et al. 2000; Kataoka et al. 2001; Chatterjee et al. 2008). 
Hence, the variation amplitudes reported in subsection 3.4 
are considered to be those of the long-term variation com- 
ponent. In subsection 3.2.3, we reported that the color 
of the variation component was constant in most objects. 
Since this constant-color feature is a result from the whole 
data set obtained during our observation period, it is prob- 
ably a feature of the long-term variation component. 

Thus, the implications for the blazar variability ob- 
tained from our observation are summarized below: The 
short-term flares on time scales of days-weeks tended to 
exhibit spectral hysteresis; their rising phase was bluer 
than the decay phase around the flare maxima. They 
were occasionally superimposed on the long-term varia- 
tions on time scales of months-years. The color of the 
long-term component was apparently constant. In addi- 
tion to those short- and long-term synchrotron compo- 
nents, a stationary component was probably embedded. 
The bluer-when-brighter trend is commonly observed be- 
cause this stationary component had a redder color than 
the variation components. A part of FSRQs has a blue 
and thermal emission as an underlying component. 

4-2. Physical Implication from the Bluer-When-Brighter 
Trend 

Several models have been proposed to explain the bluer- 
when-brighter trend in blazars. When we consider two 
distinct synchrotron components, the bluer-when-brightcr 
trend could be explained if a flare component has a higher 
fpeak than an underlying component. Even when we con- 
sider a single synchrotron component, it could be ex- 
plained by a shift of its f pea k to a higher region. In either 
case, high z/ pea k can be obtained by the injection of en- 
ergy into emitting regions, and thereby, the increase of 
the number of high energy electrons. The most plausible 
scenario for the origin of the energy injection is internal 
shocks in rclativistic shells (e.g. Spada ct al. 2001; Zhang 
et al. 2002). Besides the energy injection scenario, the 
blucr-when-brighter trend can be explained by a change 
in the beaming factor, 6, of the emitting region. This is 
because it leads to not only an apparent increase in the 
observed flux, vf v cx <5 4 , but also a shift of the observed 
frequency, ^ b s oc 8 (e.g. Villata et al. 2004; Papadakis 
ct al. 2007; Larionov et al. 2010). 

The short-term variations on a time scale of days-weeks 
tended to exhibit the spectral hysteresis in the Fj-Fy dia- 
gram, as reported in sub-subsection 3.2.3. This hysteresis 
presumably has the same nature as the hysteresis which 
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have been seen between the flux and color, or spectral in- 
dex (Kataoka et al. 2000; Ciprini et al. 2003; Wu et al. 
2007b). Kirk et al. (1998) reproduce the hysteresis be- 
tween the flux and spectral index using their model for 
the acceleration of electrons in a shock region. The fact 
that the short-term variations tended to exhibit the spec- 
tral hysteresis supports the energy injection scenario. It is 
difficult to explain such a hysteresis pattern with a change 
in 5. 

The long-term variation components on a time scale 
of months-years apparently had a constant color. The 
bluer-when-brighter trend was observed because the un- 
derlying source was redder than the variation component, 
as mentioned in sub-subsection 3.2.3 and the last subsec- 
tion. The emission of the underlying source is presumably 
synchrotron emission, because the thermal emission from 
AGN should be much bluer than the synchrotron emis- 
sion of the long-term variation component. The red color 
of the underlying source suggests that its f pea k is lower 
than that of the long-term component. 

In subsection 3.4, we showed that a higher ^ pGa k objects 
exhibited smaller variations all in the flux, color, and PD. 
Based on the discussions in the last subsection, it suggests 
that the variation amplitude of the long-term component 
was small in a region of v < f pea k and vice versa. This is 
preferable for the scenario that the long-term variation is 
also attributed to the energy injection events. If changes 
in 5 would be the major cause, it means an eccentric sit- 
uation that S changes only occur in components having 
high fc'pcak- Otherwise, large variations would be observed 
in a region of v < z^ poa k when low f pea k components are 
amplified by an increase in S. The energy injection sce- 
nario can explain the f pea k dependence of the variation 
amplitude without such an eccentric situation. 

4-3. Implications from Polarization Variations 

The high fraction of bluer- when-brighter objects implies 
that there is a major mechanism common to blazar vari- 
ations. The most plausible explanation is the energy in- 
jection scenario, as mentioned in the last subsection. The 
internal shock is a promising candidate for the mechanism 
of the energy injection to the emitting region. If the in- 
ternal shock results in an aligned magnetic field in the 
emitting region, a universal aspect of polarization varia- 
tions can be expected. 

In contrast to the strong correlation between the flux 
and color, the correlation between the flux and PD was 
weak; 30 % and 12 % of the sample exhibited significant 
positive and negative correlations, respectively. On the 
other hand, these fractions of objects are too high to con- 
clude that the polarization variations are completely ran- 
dom. It is possible that there is no universal law of po- 
larization associated with blazar variations. If this is the 
case, our observation indicates that flares could occur in- 
dependent of the structure of the magnetic field in the 
flaring source. 

Alternatively, a universal law could be hidden in the ob- 
served polarization, although it was present. Polarization 
can be described with two parameters, for example, PD 



and PA. It can also be treated like a two-dimension vector 
with Stokes QU parameters. This property of polarization 
is an important difference from that of color. The weak 
correlation of the flux and PD may be caused by the pres- 
ence of two or more polarization components. Even if a 
flare has specific PD and PA, they would disappear in the 
observed PD and PA variations because of the superpo- 
sition of other flares having different PD and PA. Based 
on this picture, we have developed a method to separate 
the observed polarization into two components: A short- 
term variation component having a positive correlation 
between the total and polarized fluxes, and a long-term 
one (Uemura et al. 2010). Applying this method to our 
blazar data, we found that the behavior of polarization in 
OJ 287 and S2 0109+224 can be explained with this two- 
component picture. The two-component picture is also 
supported by the color behavior as described in subsec- 
tions 4.1 and 4.2. 

There are three well-observed blazars in our sample 
which showed a significant negative correlation between 
the flux and PD: OJ 49, ON 231, and 3C 66A. In ad- 
dition, temporary negative correlations were observed in 
BL Lac and OJ 287, as mentioned in subsection 3.3.2. It is 
noteworthy that the observed QU of all those five objects 
do not center at the origins oiQU , but concentrate around 
regions significantly apart from the origins. Such distri- 
butions are also seen in their (Q/I,U/I) plane. These 
distributions strongly suggest the presence of long-term 
polarization components. The negative correlation means 
that the flux variations were associated with the system- 
atic changes in polarization. If the long-term component 
is subtracted, the variations might change to those show- 
ing positive correlation. Hence, the negative correlation 
can also be explained by the two-component scenario. 

5. Summary 

We performed photopolarimetric monitoring of 42 
blazars in the optical and near-infrared band using the 
Kanata telescope. Our findings are summarized below: 

• The bluer-when-brighter trend was observed in 28 
blazars, which correspond to 88 % of our well- 
observed sample. Our observation unambiguously 
confirmed that it is a universal aspect in blazars. 

• The redder-when-brighter trend was observed in 
three objects (3C 454.3, PKS 1510-089, and 
PG 1553+113) when they were faint. Even in those 
three objects, the bluer-when-brighter trends can be 
seen when they were bright. 

• The bluer- when-brighter trend was, in general, gen- 
erated by a variation source apparently having a 
constant color and an underlying red source. 

• Prominent short-term flares tended to exhibit a 
spectral hysteresis; their rising phase was bluer than 
their decay phase around the flare maxima. 

• Significant correlations of the flux and the polariza- 
tion degree were observed only in 10 blazars, which 
correspond to 30 % of our well-observed sample. 
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The correlation was weak compared with that be- 
tween the flux and color. 

• The rotations of polarization were observed in 
PKS 1510-089, 3C 454.3, and PKS 1749+096, and 
possibly in S5 0716+714. We demonstrated that 
false events of polarization rotations could be de- 
tected in the case that multiple polarization compo- 
nents are present. 

• We found that the variation amplitudes were smaller 
in low ^p Ca k objects all in the flux, color, and polar- 
ization degree. 

Based on those results, we propose that there are several 
distinct variation sources on different time-scales, colors, 
and polarizations in blazars. Both short- and long-term 
variation components are probably attributed to the en- 
ergy injection into emitting sources, for example, with the 
internal shock in relativistic shells. 

This work was partly supported by a Grand-in-Aid from 
the Ministry of Education, Culture, Sports, Science, and 
Technology of Japan (22540252 and 20340044). 

Appendix 1. Temporal Variations in Flux, Color, 
and Polarization 

Figures 21-34 show temporal variations in flux, color, 
and polarization of blazars obtained by our observation. 
The left panels are temporal variations in five parame- 
ters. From top to bottom, the panels show the V~-band 
light curve, the V — J color, the polarization degree (PD) 
in percent, the polarization angle (PA) in degree, and the 
corrected PA. The last one was corrected for the 180° am- 
biguity of PA (for a detail, see sub-subsection 3.3.3). The 
upper and lower middle panels show the color-magnitude 
diagram and the PD-magnitudc diagram, respectively. 
The upper and lower right panels show the distribution 
of (Q/I,U/I) and (Q,U), respectively. 
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Fig. 23. PKS 1510-089 (top), PKS 1749+096 (middle), and OJ 287 (bottom). 
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Fig. 24. 3C 273 (top), AO 0235+164 (middle), and OJ 49 (bottom). 
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Fig. 26. PKS 0754+100 (top), BL Lac (middle), and S4 0954+65 (bottom). 
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Fig. 30. 3C 66A (top), PG 1553+113 (middle), and ON 325 (bottom). 
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Fig. 31. PKS 0422+004 (top), H 1722+119 (middle), and PKS 2155-304 (bottom) 
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